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Abstract 

In this paper we show how using a relativistic kinetic equation the ensuing expression for 
the heat flux can be casted in the form required by Classical Irreversible Thermodynamics. 
Indeed, it is linearly related to the temperature and number density gradients and not to 
the acceleration as the so called "first order in the gradients" theories propose. Since the 
specific expressions for the transport coefficients are irrelevant for our purposes, the BGK 
form of the kinetic equation is used. Moreover, from the resulting hydrodynamic equations 
it is readily seen that the equilibrium state is stable in the presence of the spontaneous 
fluctuations in the transverse hydrodynamic velocity mode of the simple relativistic fluid. 
The implications of this result are thoroughly discussed. 
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I. INTRODUCTION 



In his semina 



Eckart in 1949 [l] introduced a bold assumption which has had a profound influ- 
ence on the subject. Indeed, in order to account for dissipative effects in a non-ideal 
(viscous) fluid he introduced the heat flow present in the fluid as part of Einstein's 
matter-stress tensor T M „. Being heat a non-mechanical form of energy, this has 
raised a rather strong debate about the physical meaning of such an assumption 
Moreover, his formalism often referred to as a first order theory jf], Q] has 
given rise to additional controversies. In a classical work Hiscock and Lindblom [8| 
showed that in the resulting linearized version of the hydrodynamic equations, the 
fluctuations in the transverse mode of the hydrodynamic velocity grow exponentially 
with time. This is a consequence of the fact that Eckart, in view of his assumption, 
was rather forced to introduce a constitutive equation for the heat flow Jj^j given by 
the equation 

J[Q] = "« (*°% - (1) 

where k is the heat conductivity, T is the local temperature, hp = 5^ + -^u a up 
is the spatial projector, u a is the hydrodynamic velocity and u a its proper time 
derivative. Also c is the speed of light and u a u a = — c 2 . Notice further that Eq. (DQ) 
is at odds with the basic tenets of linear irreversible thermodynamics (LIT). In fact, 
this formalism requires that constitutive equations linearly relate the forces with the 
fluxes [9], [jl3]. This requirement is met by the first term in Eq. (pp) but certainly 
not by the second where neither T/c 2 nor ii a represent a thermodynamic force or a 
flux. Further, physical insight of this feature of Eckart's formalism has recently been 
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discussed elsewhere 



III. 



The purpose of this work is to show that Eq. (pp) does not necessarily follow 
from a relativistic kinetic theory treatment. Moreover, we establish an expression 



relating the heat flux and its thermodynamic forces, to first order in the gradients, 
that leads to an exponential decay in the spontaneous fluctuations for the transverse 
hydrodynamic velocity mode. To accomplish this task in Sect. II we derive the 
set of linearized hydrodynamic equations which result from a relativistic Boltzmann 
equation of the BGK type, the so-called Marie equation, and in Sect. Ill we discuss 
the physical implications that may be drawn from such set. 



II. LINEARIZED RELATIVISTIC HYDRODYNAMICS 

To obtain the above mentioned set of linear relativistic equations for a simple 
viscous fluid we start by considering the special relativistic Boltzmann equation in 
the absence of external forces: 

v a l« = J (/, /') (2) 

The term on the right hand side of Eq. (j2j) is the collision operator which in general 
is a bilinear form containing the products of the distribution function / evaluated 
before and after a collision takes place. When one introduces the probability for 
the occurrence of each one of these binary collisions, namely the cross section and 
integrates over all possible velocities for the particles, one obtains the general Boltz- 
mann equation. Here, for the sake of simplicity, the kernel is modeled by means of 

nn 

the BGK approximation [12|[!3j], namely 

f _ f (o) 

J (J, f) = -— — (3) 

T 

This approximation consists in replacing all the details of such binary collisions by 
the relaxation time r, which can be viewed then as an adjustable parameter in the 
theory. This assumption simplifies the ensuing arguments for which the full collisional 
term adds no significant changes. The molecular four velocity, v a is 



7c 



(4) 



where w is the molecular three-velocity and 7 = (l — w e we/c 2 ) is the usual 
relativistic factor. All greek indices run from 1 to 4 and the latin ones run up to 3. 
The local variables are the weighted averages of the microscopic quantities. In this 
context, the number density reads 

n = j frdv*, (5) 
and the thermodynamic average of a dynamical quantity if) is defined through: 



n 



jipfdv* 



(6) 



with dv* 
equation 
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121 ] . Now, standard kinetic theory leads to the energy balance 
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— (n (mc 2 7)) + (n (mc^w/)) , £ = 
In this expression one can easily identify the heat flux as: 



J[Q] = n ( m c 2 v £ ) = mc 2 / v i f r )dv* 



(7) 



(8) 



which vanishes if the average is calculated using the equilibrium distribution function. 
At the Navier-Stokes level, following a Chapman-Enskog expansion, the distribution 
function can be written as 

/ = / (0) (1 + 0) (9) 

where is the first order correction in the Knudsen parameter, a weighted measure 
of the gradients in the system. For particles of rest mass m, relativistic parameter 
z = ^2 and in the non-degenerate case, the equilibrium function reads fl3. 
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(10) 



in which k is Boltzmann's constant, and K, n is the modified Bessel function of the 
77.-th kind. Notice also that, if the fluid is at rest, vPvp = —7c 2 . 

The function is obtained by a simple procedure. Direct substitution of Eq. 
in Eq. {2} leads to 



« f (o) a fdf W S/TO ,0/«> , 



where the second equality is due to the fact that the distribution function is a time 
dependent functional of the local variables n, T and u u . The derivatives in Eq. (TTT1) 
are given by 



#f(0) y(o) 
<9n n 



(12) 



1 + M|L_ MIL^ (13) 



^ 2 2^/C 2 (i) 2^/C 2 (J) y T 

and 



d_f^ = v l _ /(0) 
9m' 3 zc 2 



(14) 



Notice that a hydrodynamic acceleration contribution appears in the last term of 
Eq. ([TO for a = 4. If this term is left untouched and the heat flow is calculated 



with Eq. ([8]) one is led to the difficulty inherent in Eq. 



U). On the other hand, 



3-Q 



this acceleration 



according to the prescription of Hilbert-Enskog's method 
term is eliminated using Euler's equations, which precisely relate the acceleration to 
the pressure gradient. Moreover, this method establishes that by choosing the local 
variables, n, T and u v defined only through /(°), one guarantees the uniqueness of 
the solution given in Eq. Consistently with the fact that we have chosen n, 

T and u v as the state variables describing the local thermodynamic states of the 
fluid, Vp has to be expressed in terms of them. This is accomplished using the 
local equilibrium assumption which allows expressing Vp in terms of VT and Vn in 
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complete agreement with the tenets of classical irreversible thermodynamics namely, 

Ma vT +(IX v " (i6) 

Finally, since for a simple inviscid fluid whose dynamics is described by Eq. (j2J), the 
equation of state p = nkT holds true locally a straightforward calculation leads 
to the following form for the heat flux 

J IQ] = ~ l t^t ~ L n— (16) 

Equation ( TT6l ) is of the canonical form required by LIT where the thermodynamic 
forces are the temperature and the density gradients, respectively. The transport 
coefficients are given by 



-^L mc r J I 1 \ f f (0)^2 U,2 ^3/2 



/<V ( 7 2 -i) 3/2 rf 7 



and 



r 47r 7 / /" 

L n = r— mc <^ / 



/<») 7 (t> - I) 3 ' 2 d 7 - / /<V (t 2 - if 2 d-, 



(18) 

where e is the internal energy per particle. The integrals in Eqs. (IPTI ) and (fT8l) can 
be calculated in the comoving frame. Thus, the transport coefficients can be written 

as 

L T = nmc 4 z 2 Tf T (z) L n = nmc 4 z 2 Tf n (z) (19) 

where the functions fr and /„ are given by 



z 



+ 2,/C 2 (i) + 2,/C 2 (i)h/C 2 (I) {20) 



As has been shown elsewhere [15| . in the non-relativistic limit, fr{ z ) — > 5/2 and 
/ n (z) — > 0. Thus, the classical value Lt = f nk ^ 2 T is recovered for the thermal 
conductivity while the second term in Eq. (fl~6l) vanishes rendering it a completely 
relativistic effect 

The derivation of the relativistic hydrodynamic equations using Eq. (1T61 ) as the 
constitutive equation for the heat flux is a standard one and has been given in detail 
in various papers using Fourier's equation [si, [ll]. The additional term in Vn poses 
absolutely no problem at all, so we simply state the final form for their linearized 
version namely, 

5h + n 59 = (22) 
4 (rio^o + Po) 5ix v + —5n >v + ——5T V 

-(56, u - 277 (*r v % - ^5T V - ^6n, v = (23) 



nCJT + \^\ 69 - (L T 5T' k + L n 5n' k ), k = 



(24) 



where 9 = u? a , ( and r] are the bulk and shear viscosity coefficients respectively, t£ is 
the traceless symmetric part of the velocity gradient tensor, the isothermal com- 
pressibility, (3 the thermal expansion coefficient and C n the specific heat at constant 
particle density. Nought subscripts denote equilibrium quantities and 59, 5T and 
5u v denote the spontaneous fluctuations of the state variables (59 is related to 5n 
via Eq. (1221) ) around the equilibrium state of the fluid. 

We remind the reader that this set of equations is the relativistic analog of the one 
obtained from classical irreversible thermodynamics to test the validity of the linear 
regression of fluctuations hypothesis introduced by Onsager over seventyfive years 
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so HE? 



ago Il7|,ll8fl as a basic ingredient required to derive his famous reciprocity relations. 
The fact that this assumption has been experimentally tested for a variety of systems 
{19! ] calls immediately for its relevance in the relativistic regime. To illustrate this 
point let us consider the fluctuations associated to the transverse mode of the velocity 
field u v . This mode is easily uncoupled from the rest of the fluctuations by taking 
the curl of Eq. (|23|) . Since all the gradient terms disappear we are led to the result 
that 

\ (n s + po) 5U V - 2r,V 2 5U u = (25) 



where 5Uj = A k u% is the transverse mode. Taking the Fourier transform of Eq. ( 1251) 
we get that 

SU V (k,t)= 5U U ( k, O] Exp ( 27]k ° ) (26) 

V / \ J \ n Q e +p J 



where 8U V yk, tj denotes the Fourier transform of 5U U . Since the exponential's argu- 
ment is always negative, the fluctuations die out with a characteristic time ^^gr^ ■ 
In the classical limit, since ! ^f a — ► p an d § is negligible, the decay time reduces to 
2^p- in agreement with classical hydrodynamics. This result implies that the sponta- 
neous fluctuations of such mode obey the linear regression hypothesis. Notice that, 
since the last two terms in Eq. ( |23l ) are precisely Jjgj, this exponential decay is due 
to the fact that the curl of this vector vanishes. Thus this result is independent of 
the equation of state (see Eq. (fT5ll ). 

From the remaining equations, Eq. (!2~2~I) . the divergence of Eq. (l23l) and Eq. (1241) . 
one can get a set of coupled equations for 5n and 5T whose solution is related to the 
well known Rayleigh-Brillouin spectrum. A full discussion of this feature is too long 
and technical to be presented here and will be published elsewhere H|. 
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III. DISCUSSION AND CONCLUDING REMARKS 



Examination of the results obtained in the previous section manifestly indicate 
that Eq. ( fT6l ) is the main result of this paper: it satisfies all the requirements 
imposed by classical irreversible thermodynamics and leads to a set of linearized 
hydrodynamic equations which satisfy Onsager's linear regression of fluctuations hy- 
pothesis. The crucial step was to eliminate the acceleration term appearing in Eq. 
dTTT) in terms of the pressure gradient through Euler's equations. This is a point 
;hat was overlooked by Eckart himself and later by some of his followers and critics 



IB 



Thus the term ii v , which is not a force, has been replaced by a combination of 
gradients of the state variables. This substitution is necessary in order to arrive at a 
closed set of evolution equations by further using the constitutive equations, which 
are written in terms of such gradients. This is indeed the fundamental purpose of 
the constitutive equations. If u u is considered a state variable the description of 
the system is overdetermined. This is precisely the main source of the objections 
raised against Eckart's theory which materialized in the appearance of the so called 



generic instabilities 



211 ] and presumably in giving rise to transport equations of 



the parabolic type. This last aspect however deserves a much closer attention and 
will not be addressed here. In fact, the curious reader may easily verify that if u a 
in Eq. (DQ) is expressed in terms of Vp through Euler's equation all the difficulties 
mentioned above will be eliminated. 

In kinetic theory things are quite different. Although, as mentioned in the text 
an acceleration term in Eq. (fTTT) is present, it is automatically disposed off when one 
appeals to the Chapman-Enskog method for solving Boltzmann's equation. Indeed, 
when Eq. ([9]) is substituted in Eq. j2]), the term to order zero in the gradients 
yields Euler's equations. Successively, the equation for the correction of order n in 
the gradients will have a solution if and only if the hydrodynamic equations of order 



9 



n — 1 in the gradients hold true. This is the spirit of the original method of Hilbert to 
solve Boltzmann's equation and this is why the drift term in the n th order equation 
which is already of order one in the gradients has to be computed with the solution of 
order n — 1. This explains why the use of Euler's equations is demanded to calculate 
(ft in Eq. f TTTTl . We should also emphasize that these results have been shown to hold 
for a relativistic Boltzmann equation with its full collision term [22| and further, from 
phenomenological arguments, to be valid for any arbitrary relativistic fluid [111 ]. 

Moreover, to fully grasp the relevance of this calculation we can compare Eq. 
f[25T > with its counterpart obtained in Ref. [8|, in which the coupling of heat with 
acceleration is retained. In that case, the transverse velocity mode fluctuations satisfy 
the equation: 

-±8U V - -» (n e + po) 8U V + 2 V W 2 5U U = (27) 
c 4 cr 

After taking the Fourier-Laplace transform in this equation, we now obtain two 
roots governing the time behavior of the fluctuations. The calculation in a specific 
example, water at 300 K at one atmosphere leads to an extremely short characteristic 
time in the growing mode. The conclusion included in that paper states that first 
order theories should be discarded in favor of more complicated formalisms. Here we 
have shown that the apparent trivial replacement of acceleration with the pressure 
gradient while solving Boltzmann's equation, with the Chapman-Enskog method, 
leads to a nice exponentially decaying transverse velocity mode. Thus, the so-called 
generic instabilities are indeed absent in first order relativistic hydrodynamics. In 
addition to this application, we can assert that the calculation of the Rayleigh- 
Brillouin spectrum using the acceleration term leads to the unphysical result that 



such spectrum is inexistent. This has been fully discussed in Ref. 



20]. 



There is a last comment which is pertinent in view of the results here presented. 
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As we said in the introduction, Eckart's formalism as well as all others which con- 
tain an entropy balance equation with an entropy flux is given by JLJT, have been 
referred to as "first order theories". Due to the so-called generic instabilities and the 
apparent violation of the principle of antecedence they where substituted by second 



order theories 



22l | which postulate that both the entropy density as well as 
the entropy flux should be functions of the heat flow and the viscous tensor. The 
point here is to remind the reader that the first order theories, if the above defini- 
tion is kept, together with the ordinary linear force-flux relations lead, for a simple 
viscous fluid in the non-relativistic case, to the Navier-Stokes-Fourier equations of 
hydrodynamics. These equations in their full expression consist of a set of five non- 
linear coupled equations (first order in time derivatives and second order in spatial 
derivatives). Their linearized version is in total agreement with Onsager's regression 
assumption and only upon some drastic simplifications lead to transport equations 
of the parabolic type, such as Fourier's heat conduction equation. To assert that 
in the relativistic domain such theories are incomplete and/or incorrect is equiva- 
lent to asserting that there is no relativistic extension of the NSF equations. What 
we have proved in this paper is that this requires further study. What is the full 
and correct version of these equations may still be today a problem which has not 
been completely solved but already some attempts have been clearly discussed in the 
literature fl. 

We believe that these considerations are important due to the enormous effort 
that has been made in the past few years to describe the hydrodynamics of the hot 



dense matter produced in heavy ion collisions. Some workers in this fie 



the model proposed by second order theories or modified versions of it 



require the use of adjustable parameters 



d are using 



25fl-[3jj which 



As an example of another alternative 



to these theories there is Landau's hydrodynamics for hot dense matter based on 



Euler's equations. In an elegant paper by C. Y. Wong 
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331 ] . this theory has been 



recently shown to provide a very good explanation for the collective dynamics of 
the fluid that is produced in the relativistic collisions of heavy ions. If one would 
consider adding dissipative effects it is clear that the NSF equations would be the 
best candidate. Thus we conclude paraphrasing S. Weinberg's words "it is worth our 
while to develop the outlines of the general theory of a relativistic imperfect fluid" 

a 
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